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The mean spectral reflectance over large space/time scales is required 
for many climate studies (e.g., the attribution of climate change signal and 
the perturbation analysis of a climate system); CLARREO has proposed to use 
it as the climate benchmark spectrum. 

This mean spectral reflectance can be derived from satellite 
measurements or from model simulation, but measured RS spectrum over 
globe is limited. 

Spaceborne sensors have measured the atmospheric and surface 
properties over globe for decades and have accumulated a large volume of 
instantaneous measurement data that could be used to simulate the spectral 
radiation over long time and globe. However, it is computationally too 
expensive to simulate the mean spectral radiance over large climate domains 
by explicit RT computations at the satellite footprint scales. 

A simplifying approach for efficient RT computation has to be 
developed to use the instantaneous satellite data. 
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Inputs: MODIS 
clouds, aerosol, 
GOES5 water vapor 
in CERES SSF; 
SMOBA ozone; 
SeaWiFS Chi. 
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An example of model-observation comparison of spectral reflectance. This 
example is for one zone (30S-35S) in the January of 2006. (ocean only) 

Input data are averaged in each month: 5x5 deg grid ^ 5 deg latitude zones. 
Model computation in each zone: land, ocean ^ clear, cloudy ^ ice and water. 
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An example of the probability distribution function (PDF) of the instantaneous 
cloud i in one latitude zone for months of April. A different color in each panel 
represents a different year. 


Based on the cloud x PDF, we partition the cloudy footprints in each zone 
into a number of bins to account for the wide variation of cloud properties in 
different satellite footprints. 

The PDF approach essentially organizes or groups the large number of 
satellite instantaneous measurements according to the x value, so that those 
footprints with similar x and other atmospheric/surface properties are treated 
once, rather than footprint by footprint, in the RT modeling. Therefore, the 
computation time is reduced significantly compared with the computation 
footprint by footprint. 

Results for every x bin are then summed and averaged to obtain the 
mean radiance with the weight in each bin equal to the integral of the PDF in 
the bin. 


Jin, Z. et al. (2013), An efficient and effective method to simulate the earth spectral 

reflectance over large temporal and spatial scales, Geophys. Res. Lett., 40, 374-379, doi 
10. 1002/grl. 50116. 
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Comparison of observed and calculated monthly and globally averaged spectral 
nadir reflectance (ocean only). 


Latitude 


The PDF algorithm is tested in the five large latitude regions. 


Five latitude regions 
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The RMS model discrepancy in the monthly mean reflectance in the five latitude 
regions (based on the results in the 72 months from 2004 to 2009). Each color 
bar represents a RMS error for a different modeling approach. 


Application to Solar Fingerprinting Test 


The climate change fingerprinting is to attribute the large spatiotemporal 
averaged spectral variation to different climate variables: 
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AR Reflectance change signal (spectral difference between two climate states) 
K Kernel matrix (fingerprints), has to be generated from RT modeling 
AX Climate variable changes to be retrieved 
e Errors or residuals that cannot be explained by fingerprints 


Solution by least squares estimation (LSE): 


AX = (K T K)~ l K T AR 



Solution by optimal detection: 


A X = (K t E~ 1 K)~'K t E~'AR 


(2b) 




The formulation for fingerprinting retrieval here is similar to those used in 
many conventional retrieval methods applied to instantaneous satellite data. 
However, the fingerprinting retrieval differs from the traditional remote sensing 
retrieval in that 

1) It uses the average-then-retrieval approach and thus is associated with the 
averaged quantities over large spatiotemporal scales instead of the local or 
instantaneous values. 

2) It is for the difference (AR and AX) between two climate states instead of the 
absolute value (R and X). 

The linear expression above implies that AR and AX have to small to maintain sufficient 
linearity. 


We use model-simulated data to test the concept of fingerprinting 
attribution. Model input parameters: 

S 10 years (2001-2010) of CERES SSF data from NASA Terra satellite ; 

S aerosol and cloud properties (optical depth, particle size, phase, and 
height) retrieved from MODIS; 

S column water vapor and surface wind data from GEOS5-MERRA 
reanalysis; 

S ozone data from SMOBA; 

S ocean chlorophyll concentration from SeaWiFS. 

Input these parameters to the combined COART-MODTRAN model, we 
generated the spectral kernels and a time series (10 years) of monthly mean 
reflectance spectra (320-2300nm, 4nm resolution) to retrieve the interannual 
changes in the 11 relevant climate parameters: 

PW, AOD, 0 3 

and 

T, Fc, Ht, Re for water and ice clouds, respectively . 
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An example of the model-simulated monthly and global mean reflectance 
anomalies (AR) in the four months. In each panel, each color is for a different year. 


AR is indeed small: typically less than ±3% of the mean reflectance 
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An example of solar spectral reflectance kernels 


Reflectance Kernel (1/c 
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Comparison between the 
fingerprinting retrieval and the 
observational truth for the 
monthly global mean anomaly. 

Each panel is for a different 
climate parameter (11 total). 

The nonlinear error in kernels is 
not considered in the retrieval. 


Retrieval 
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Same retrieval as above, but 
the nonlinear error in kernels is 
considered. 


Two nonlinear error sources: 


1). nonlinear radiative response : 


f 


x 0 +Ax 


K(x)dx * K(x 0 )Ax 


2). radiative interactions: 


AR*Z i AR j =Z i KAX i 
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Without nonlinear errors considered 
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With nonlinear errors considered 
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Without nonlinear errors considered 
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With nonlinear errors considered 
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RMS comparison in all regions 
fro all 11 parameters. 

Each panel represents a 
different region. 

The RMS is normalized by the 
variation a for each parameter. 


■ Without nonlinear error considered (LSE) 

| With nonlinear error considered (optimal detection) 


Conclusion 


S The cloud PDF approach provides a simple, fast, and effective way to simulate the 
mean spectral reflectance over large time/space scales with a large volume of 
satellite data. 

S Using ten years of satellite data and the PDF method, we tested the concept of 
using fingerprinting approach for climate change detection/attribution. 

S Comparing the fingerprinting retrieval to the observational truth, the RMS 
differences are less than 2 a of the variance for all variables in all regions. A large 
error usually corresponds to the variables with large nonlinear radiative response, 
such as cloud x and ice particle size. 

S Using the optimal detection method to take into account the nonlinear radiative 
error in the kernels, the retrieval accuracy is significantly higher, so that the RMS 
errors are reduced to less than la of the variance, indicating the profound impact 
of the nonlinear error on fingerprinting retrieval. 

S Ozone retrieval has the highest accuracy among all variables in all cases. 

S The test results demonstrate that the concept of climate change fingerprinting 
based on the reflected solar benchmark spectra is viable. 


